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Abstract 

We present a new and simple bound for the exponential decay of 
second order systems using the spectral shift. This result is applied 
to finite matrices as well as to partial differential equations of Mathe- 
matical Physics. The type of the generated semigroup is shown to be 
bounded by the upper real part of the numerical range of the underlying 
quadratic operator pencil. 



1 Introduction and main estimate 

In this note we consider the abstract second order system 

Mx + Cx + Kx = 



(1) 



where M, C, K are selfadjoint operators in a Hilbert space X with the scalar 
product x*y linear in the second and antilinear in the first variable.^ For 
simplicity, we assume that the operators M, K, C are positive definite and 
that M, C are bounded. The phase space formulation of reads 



y = Ay, A 







K'/^x 



(2) 
(3) 



with the solution 



y 



Thus, the square of the norm equals twice the total energy of the system: 

|2 II 7>'1/2„||2 



\y\ 



\K^/^xr + x*Mx. 



*Fernuniversitat Hagen, Lehrgebiet Mathematische Physik, Postfach 940, D-58084 Ha- 
gen, Germany, e-mail: kresimir.veselic@fernuni-hagen.de. 

^AU other conventions and notations will be taken from [S]. 
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The operator A is readily seen to be maximal dissipative on 



and thus the semigroup e"^* is contractive i.e. the energy of the system is a 
non-increasing function of t? The resolvent of A is immediately seen to be 
given by 



(4) 



{A-\I)-^ = 
with 

K{\) = + \C + K (5) 

at least for those A for which K{X) remains positive definite. This formula is 
rigorous, if M, C, K are all bounded, otherwise the terms /'^ K {X)~^ /"^ 
and K^^'^K{\)~^K^/'^ have to be replaced by their closures; the latter are 
obviously everywhere defined and bounded. 

Most existing works on the exponential decay estimate the infimum of all 
(3 for which 

lle^i <C^e^*. (6) 

The infimum value ojq{A) of all possible /3 in © (the type of the semigroup 
in the terminology of |i5j) is often equal to the maximal real part of the spec- 
trum of A. In looking for this infimum usually little attention is paid to the 
constant Cp which may tend to infinity as (3 approaches the infimum (see e.g. 

[2]; in] )• Since this constant plays a key role in controlling the finite-time 
behaviour of the system we are interested in a bound in which both j3 and 
C/3 are tried to be made simply computable from the coefficients M, C, K. 
In fact, we obtain a family of estimates (jHl) for any j3 from the interval (7,0], 
where 



-x*Cx + Jix*CxY - Ax*Mx\\K'^l'^x¥ 
7= sup 3? ^— — -— (7) 

while Cp is an expression with similar ingredients. The set 

W{K) = {A G C; x*K{X)x = 0, for some unit x} 

is called the numerical range of the matrix pencil K{\) (cf. e.g [Hj). Our result 
implies, in fact, 

uJoiA) < 7 = snp ^W{K). 

^The reader primarily interested in finite dimensional applications may skip operator- 
theoretical details and just take M, C, K as matrices. In this context 'positive' means 
'positive semidefinite' and the expression ||if^/^a;|p below may always be read as x*Kx. 
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The bound for ujq^A) , obtained in [2] reads in our notations 



, „ X Cx 1 
7b = max < — mi — — , = 

> X 2x*Mx x-Cx I o x*Kx 

The values 7 and 7^, are not easy to compare in general. For any underdamped 
system i.e. whenever the expression under the square root in ((7j) is uniformly 
negative, we obviously have 7 < 7^. For further comparisons see Sect. 2 below. 

An estimate for Cp was obtained in [3', Vol. I, Ch. 6. for the wave equation 
with distributed viscous damping. The bound obtained there is much less 
explicit than ours; this actually made impossible any comparison of the two. 

As a by-product, we prove that uJq{A) is equal to the supremum of the real 
part of the spectrum of A for "partly overdamped systems" i.e. for those for 
which 27M + C is positive definite. Another case in which this equality is 
shown to hold are the so-called modally damped systems. 

Our main tool will be the 'spectral shift substitution' 

X = e^'z 

which gives rise to a new phase space representation, equivalent to the previous 
one. This yields a fairly simple total-energy decay estimate. Applications 
are made to both finite matrices and differential operators of Mathematical 
Physics. The obtained estimate is shown to be (asymptotically) attainable 
and therefore in some weak sense optimal. On the other hand, our estimate is 
void, if damping has a nontrivial nulspace. While there are such systems with 
no exponential decay at all, there are still many relevant cases whose decay is 
not covered by our theory. In Sec. 2 we provide illustrating examples. 
The substitution x = e^^z yields 

X = jje^^z + e^*i 

X = /i^e'^*^ + 2/xe^*i + e>'^z 

and ^ reads 

Mz + C(/i)i + K{fx)z = (8) 

with 

C{fi) = 2/iM + C. (9) 

As long as C{fi),K{fi) stay positive definite this is equivalent to the phase 
space representation 



w = Aw, A 



ir(/i)V2M-l/2 



(10) 
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w 



with the solution 



'Ill 



We now connect these two representations. We have 



yi 



Thus, 



Ky^K{ii)-y^ 
/iMV2K(;^)-i/2 J 

This, together with the evolution equations for w gives 

C{^^)A={A-^^I)C{^^). 

This yields the decay estimate 

||e^*|| < ||£(/i)||||£(/i)-ie^*. 



(12) 
(13) 

(14) 
(15) 



The foregoing formal calculation is rigorous for finite matrices. Our general 
operator setting requires additional justifications of these steps. Also, we need 
a more explicit bound on the condition number which appears on the right 
hand side of (fTHjl . These two issues are the subject of the following. We set 



x*Cx 
2x*Mx 



± jD(x) 



(16) 



th^ 



wi 



and 



, / x*Cx \2 \\K^/^xP 

Dix) = — 

^ ' \2x*MxJ 2x*Mx 



7 = sup 3?p+(x) 

a;gD(ifl/2) 



(17) 



Proposition 1 Let M,C,K be as defined above. Then (i) 

„ x*Cx 
7 > — mf 



^ 2x*Mx' 



(18) 



moreover, 7 is the infimum of all fi for which both K{fi) and 2/iM + C are 
positive definite. 

■^By convention, the square root of a negative number will have the positive imaginary 
part. 
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Proof. The relation (|18|) is obvious. To prove the second assertion note that 
for any real fi the selfadjoint operator K{fi) is generated by the form 

hf,{x, y) = iJ,^x*My + fix*Cy + {K^/^x)* K^''^y 

which is symmetric, bounded from below and closed on 'D{K^/'^). Take any 
/i > 7. Then 

h^{x, x) = — — p+{x)) > x*x{fj. — 'Sip+{x)Y > x*x{jj, — 7)^ 

(note that p_ (a;) < p+{x) whenever D{x) > 0). Thus, K{fi) is positive definite. 
By fITHj) 2fj,M + C is positive definite also. 

Conversely, suppose that 27M + C is positive definite; there is a sequence 
Xn of unit vectors such that 

3^P+(^n) = -T^i^ + 3fjy^(^0-.7, ^-00. 
Since 27M + C is positive definite it follows 

lim inf ^JD(xn) > 0. 

n^oo V ^ ' 

Thus, for n sufficiently large D{xn) > 0, hence 

P+{Xn) = "^p+iXn) ^7, 00 

and (-^'(7) is positive semidefinite) 

h^,{Xn, Xn) = WKi^jy/^Xnf = {l - P~{Xn)){l " P+i^n)) ^ 

since p_{xn) is obviously bounded. Thus, -^'(7) has not a bounded inverse, 
and is not positive definite. Q.E.D. 

Theorem 1 Let M, C, K be as defined above. Assume that 7 < /i < holds. 

Then (i) the operator C{fi) together with its inverse is bounded and everywhere 
defined, (ii) Both operators leave the subspace T>{A) invariant and holds. 
(Hi) The estimate ^ holds with 

/3 = /i, Cp = \\cmmi^r% (19) 

(iv) The type of the semigroup is bounded by 

< 7. (20) 



5 



Proof. Since 7 is negative then for —7 < fi < both K{^) and C{^) remain 
positive definite. The boundedness of M,C imphes (cf jS]) 

V{K{i2y/^) =V{K^/^), and V{K{^j)) = V{K), > /x > 7, 

so both (fTT?jl and its inverse 



(21) 



are everywhere defined and bounded (a straightforward calculation shows that 
the one is the inverse of the other). This also shows that P(v4) = 'P(A). 
Take a vector from V{A), that is, x G V{K^/'^), y G M^/2-p(/^i/2) ^nd set 



Here x = K^/^K{fi)-^/'^x is from V{K^/^) because K{fx)-^/^x G V{K{fi) = 
V{K). Also V = fiM^^'^K{fi)-'^/^x + yis from M^/'^V{K^/'^) because y is such 
and K{n)~^^'^x G V{K). Thus, i2(/i) leaves ^^(^4) invariant. The proof for 
C{n)~^ is similar. Now the relation p4|) can be directly verified on any vector 
from T^IA). The estimate (|T^ is now obvious. Q.E.D. 



X 




u 






V 



Corollary 1 // 



then 



7 > — inf 



v.. X 



7 = uJo{A) = sni)a{A). 



Proof. Under our assumption, as in the proof of Proposition ^ it follows that 
K{'j) is not boundedly invertible while K{X) is positive definite for any A > 7. 
Now, the 2, 2-block in the resolvent matrix in (0} gets unbounded for A = 7. 
Thus, all A with 3?A > 7 belong to the resolvent set of A whereas 7 belongs 
to its spectrum, hence 7 = sup3?cr(^). From Theorem [T] it follows uJo{A) < 7 
while ujq{A) > supcr(^) is a general fact. So, the assertion follows. Q.E.D. 

The systems covered by the corollary above may be called 'partially over- 
damped' in a sense similar to that introduced in jl]. 

We next give some bounds for the condition number appearing in 
Note that both C{fi) and £(/i)~^ are of the type 



A 
B I 



(22) 
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Lemma 1 Let A,B be bounded operators and Then \\C\\ from is bounded 
by any of the numbers 



'1 + \\A*A + B*B\\, max{p||,l+ ||5||}, max{||A|| + ||5||, 1}. (23) 
Proof. We have 



A*A + B*B B* 
B I 



CC* 



AA* AB* 
BA* BB* + I 



If B = the assertion is trivial. If B ^ the number a = = = 
||£*£|| > 1 belongs to the spectrum of C*C and there exist sequances Xn, yn 
with ||x„f + \\ynf = 1 and 



{A*A + B*B)xn + B*yn-aXn 
Bxn + yn-ayn 



for n ^ 0. Hence 



B* B 

A* A + B*B + a]xn^O 

a — 1 



(24) 
(25) 
(26) 



(27) 



By (j2ni) the sequence x„ does not converge to and without loss of generality 
we may assume that in ll^^nll = 1 holds. Then 



Hence 



- a{{A*A + B*B + /)x„, x„) + {A*Axn, x^) 0. 



a < sup 



11x11=1 



{{A* A + B*B + I)x, x) + J {A* A + B*B + I)x, x)^ - A{x, x){A*Ax, x) 



< sup {{A*A + B*B + I)x,x) 
lkll=i 



The other two bounds in (PHj) are obvious. Q.E.D. 

Remark 1 (i) Ours is, in fact a family of bounds depending on the parameter 
from the interval (7, 0]. An optimal bound would be obtained as the infimum 
over all of them. 

(ii) By continuity, the bound (fT^ remains valid even for /i = 7, if -^'(7) is 
positive definite. In this case we have ujq^A) < 7 i.e. or bound is not optimal. 
Indeed, for some fJ^ < 'j the operator K{fj,) will still be positive definite and A 
will still generate a uniformly bounded semigroup ( 5J, Ch. IX, Th. 2.1) and 
(HI), (Hni) will still be vahd. 
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(iii) Our main estimate does not contain the norms of the operators M, C 
or their inverses and indeed it would certainly hold in much more singular 
cases, but then some additional regularity conditions would be needed e.g. 
the boundedness of the operators C^/'^K~^^'^ ^ C~^I'^M^I'^ and the like thus 
requiring more technical proofs. The same is valid, if we would admit more 
general damping operator C by merely asking it to be accretive. 

(iv) The expression (jl7p for 7 is neat but not easily computable, even in the 
case of finite matrices. As suggested by Proposition ^ a simple viable method 
to determine 7 numerically would run as follows 

• Find 7o = — inf 2x*mx finding the lowest eigenvalue of the matrix 
pencil 2AM - C. 

• If 7o ~ halt, no bound available. 

• If 7o < and -^'(70) is positive definite then 7 = 79. 

• If 7o < and -^'(70) is not positive definite then seek 7 by bisection in 
the interval (70, 0]. 



2 Examples 



Our first example is the one-dimensional system with the two-dimensional 
phase space matrix 



A 



k 
-k -d 



k>0, d>0. 



A straightforward, if a bit tedious, computation gives 



COS St + d{sm St)/ (26) k{smSt)/S 

-k{smSt)/S cos St - d{sm St)/ {2S) 



-dt/2 



(28) 



(29) 



with 



S = V4P -(^2/2. 

The above formula is valid for S both positive and negative (in the latter 
case there is a real expression by means of hyperbolic functions) whereas the 
formula for Ak"^ — d"^ = is obtained taking the limit k ^ d/2, thus reading 



I + dt/2 dt/2 
-dt/2 1 - dt/2 



-dt/2 



(30) 



Here 







^J^J?+^ld+k'^ 
, ^ 1 
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and our Theorem above yields 
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-d + Vd^ - 4fc2 



(31) 



Here the right hand side is the largest real part of the spectrum and therefore 
the inequality is, in fact, an equality i.e. our estimate for uJo{A) is sharp 
for all possible values of k, d. The bound from Pj reads (in our notations) 



76 = max 



d' 



d + 2k' 



As a straightforward calculation shows we have here 76 = 7 for d < (v^— 1)^ 
and 76 > 7 for (i > (y^ — l)k. 




Figure 1: k/d = 1 



On Figs. [11121 we display our bound (dotted line) with the true norm (solid 
line) as functions of t (our bound is obtained by taking the minimum over 
four equidistant values of /i). We see that our constant gets pessimistic for 
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Figure 2: k/d = 0.5 

small values of k/d. On the other hand, for k/d oo we may take /i = d/2, 
I and (fTTHl reads asymptotically 

||e^*|| < e''*/^ 

whereas the right hand of ()29p behaves as 

cos kt sin kt -dt/2 
— sin kt cos kt 

Thus, or bound (fTU|) is asymptotically sharp. 

Any modally damped system i.e. a system in which M, i^', C satisfy the 
relation 

CM^^K = KM^^M 

is easily seen to be unitarily equivalent to an orthogonal sum of matrices of 
type (j28j) . a bound is obtained as the maximum over all of them. The type 
ujq{A) is the maximum of the values in (jHTj) .^ In particular, ujq{A) is equal to 

As a second example take the wave equation in a bounded domain f2 C i?" 

wtt{x, t) + c{x)wt{x, t) + Aw{x, t) = 0, (32) 

*This remark is literally true if the system has discrete spectrum (e.g., if K has a compact 
inverse). Otherwise we would have direct integrals. 
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with the boundary condition 

w{x, t) = X G dil. 

The function c is supposed to have finite positive minimum and maximum. 
We will estimate e.g. the second of the bounds in ()23|) . Taking first 

A = K'/^K{fi)-'/\ B = /iMi/2^(/i)-i/2 
in X = L2{R) we have 

1 1 



\A\ 



l + inf, "'(f+^:;»)" 1 + Ai' 



where Ai is the lowest eigenvalue A of the boundary value problem 

(/i^ + fic{x))u = —\Au, u\dn = 0. 
Since under our assumptions /i^ + /ic(x) is negative definite we have 

Ai < 0, 1 + Ai>0. 

Furthermore 

2 /iVw 1 

||-d|| = sup 



1 



u* (pc(-)—A)u 
u*u 



/x^ + As' 

where A2 is the lowest eigenvalue A of the boundary value problem 

(/ic(a;) - A)n = Au, = 0. 

Taking next 

A = K{iiY'^K-^'\ B = -iiM^'^K-^'\ 
we have ||A|| < 1 and 

where A3 is the lowest eigenvalue A of the boundary value problem 

Am = —Xfiu, u\gQ = 0. 

Altogether 

= \\ci^i)\\\\c-\^i)\\ <max\ ^=L=,i + ^^^==] (1+ ^ 



VTTTi' V/u^ + A2 j V va;/' 

where Ai, A2, A3 are obtained above. Thus, our bound for is obtained from 
the extremal eigenvalues of some self adjoint elliptic boundary value problems 
involving M, C, K. If c{x) is constant then all these boundary value problems 
reduce to —Am = Am. This system is also modally damped, so the correspond- 
ing estimates are applicable here, too. 
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